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^ , Abstract 

We clarify the algebraic structure of continuous and discrete quasi-exactly solvable spectral prob- 
lems by embedding them into the framework of the quantum inverse scattering method. The quasi- 
exactly solvable hamiltonians in one dimension are identified with traces of quantum monodromy 
matrices for specific integrable systems with non-periodic boundary conditions. Applications to the 
Azbel-Hofstadter problem are outlined. 



1 Introduction 



, At present time the methods related to quantum integrability are highly developped. Originally they 

• were invented and used in the context of quantum field theory. However, the main ingredients proved to 

be purely algebraic. They can be successfully applied to quantum-mechanical (i.e., one-particle) problems 
, . ■ as well. 

rN I An example of a new application of this kind is the recent progress [Q, ||] in the famous problem of 

' Bloch electrons in magnetic field on a two-dimensional lattice Q (sometimes called the Azbel-Hofstadter 

problem). Even the one-particle problem is non-trivial. In a proper gauge it reduces to a one-dimensional 
quasiperiodic difference equation (Harper's equation is the most popular example). It has been shown 
m §, d that some of these equations admit partial exact solutions of the form typical for quantum 
integrable systems: the eigenfunctions are polynomials with the roots constrained by Bethe equations. 
These solutions give some specific states (one for each stable band), the energies being expressed through 
the roots in a simple way. In the paper this result has been generalized to a whole class of second-order 
difference operators admitting partial algebraization of the spectrum. 

The form of the result suggests to ask for a direct connection with quantum integrability. Indeed, 
such connection does exist. 

A unified approach to quantum integrable systems is most elegantly formulated in terms of the Quan- 
tum Inverse Scattering Method (QISM) created by the Leningrad School § (for a more recent review 
see 0). In the paper |^ (Appendix B) it has been shown how to embed the Azbel-Hofstadter problem 
into the QISM. The hamiltonian has been identified with a quantum transfer matrix (trace of a quantum 
monodromy matrix) for a specific integrable system with boundaries. This gives a possibility to apply 
the powerful machinery of the QISM such as functional Bethe ansatz [Q . 

This paper may be considered as an extensive comment to Appendix B of the paper j2| . We give a 
detailed construction of the quantum monodromy matrices for the general family of difference equations 
considered in 0|. In a more general context, we provide new formal grounds for studying difference or 
differential operators in one variable having a finite number of polynomial eigenfunctions. The continuum 
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limit corresponds to models with the rational i?-matrix. In this case we reproduce a class of second-order 
differential operators having the property of the partial algebraization of the spectrum. Their eigenvalue 
equations were considered in the literature some time ago [|j , ^ (the idea goes back to the papers , 
JTlf ). A systematic treatment, based on the hidden dynamical sZ(2)-symmetry, was given in [ [l^ (see 
also the reviews and references therein). A different approach was suggested in These 

equations are known as " quasi-exactly solvable" problems ("quasi" means that usually only a part of 
the spectrum can be found in a closed algebraic form). The corresponding hamiltonians are known 
to be quadratic forms in the standard generators of sl{2) (taken in a finite-dimensional representation). 
In our approach the generators of sl{2) appear as matrix elements of the universal 2x2 L-operator of 
XXX-type. To g-deform this picture, one should use XXZ-type L-operators. Their matrix elements 
are expressed through generators of Uq{sl{2)), the g-deformation of the universal envelopping algebra of 
sl{2). This construction yields difference equations [|. 

In short, we reduce the quasi-exactly solvable spectral problems mentioned above to the problem 
typical for quantum integrable systems and lattice statistical models, i.e., to diagonalization of a transfer 
matrix. Besides, we give the QISM interpretation of isospectral transformations of quasi-exactly solvable 
hamiltonians (in the continuous case) under adjoint S'_L(2)-action. 

Here is a more detailed description of the content. 

In Sect. 2 we describe universal trigonometric 2x2 L-operators depending on spectral parameter. 
There are two kinds of them: one is related to Uq{sl{2)) (it is usually used in integrable XXZ magnets 
with higher spin), another one is associated to the dual algebra, Ag{SL{2)). The matrix elements are 
expressed through the generators of Uq{sl{2)) and Aq{SL{2)) respectively. The relevant representations 
of these algebras are briefly reviewed. The rational limit of these L-operators is also discussed. The 
former turns to the universal L-operator of the isotropic XXX-type integrable magnet while the latter 
becomes a c-number 2x2 matrix independent of the spectral parameter. 

The necessary extraction from the formalism treating integrable systems with boundaries is given in 
Sect. 3. The starting point is "reflection equations" Following [Q, we recall the construction of 

quantum monodromy and transfer matrices for systems with boundaries. 

In Sect. 4 we apply this general formalism to the elementary L-operators and obtain in this way 
quasi-exactly solvable hamiltonians. Applications to the Azbel-Hofstadter problem are outlined. 

The rational (continuum) limit is treated in Sect. 5. A comparison with the representation in terms 
of Gaudin's magnet (suggested in jl^) is made. In Sect. 6 we discuss isospectral transformations of 
continuous quasi-exactly solvable hamiltonians under adjoint action of SL{2). The QISM interpretation 
of these transformations in terms of the rational limit of the L-operator related to Aq{SL{2)) is suggested. 
Sect. 7 contains conclusions and speculations on some open problems. 



2 Elementary L-operators and quantum algebras 

The standard basis of the quantum integrability is the Yang-Baxter equation ^ (YBE) with a spectral 
parameter u: 

R{u/v)Ti{u)T2{v) = T2{v)Ti{u)R{u/v) , (1) 

where Ti{u) — T{u) ® 1, T2(u) = 1 T{u). We take T{u) to be a 2x2 matrix with operator matrix 
elements and 




\ qu-q-'^u-^ ) 

is the symmetric trigonometric i?- matrix with the parameter q (or its rational degeneration). In the 
trigonometric case we use the multiplicative parametrization. 

^An attempt to apply quantum algebras for generating quasi-exactly solvable difference equations was made in [ p^ . 
However, the authors used another version of the quantum algebra (which does not allow one to construct hermitian 
hamiltonians; for details see [Q) and did not discuss the Bethe ansatz solutions. 

■^Sometimes (hi) is called the Yang-Baxter equation only if T[u) = R{u), with considering it as a cubic functional relation 
for R{u). 



q ^ 











u — u^^ 


q-q- 





q-q^^ 


u — u 












(2) 
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The equation (|lj) determines commutation relations for the elements of the quantum monodromy 
matrix T{u). The elementary solutions of the YBE (i.e., those which can not be decomposed into a 
product of simpler ones) are of particular importance for us. They are called L-operators. In lattice 
integrable models (or spin chains) an i-operator is usually associated to a lattice site. Here are two main 
examples of i-operators. 

1). L-operators associated to Uq{sl{2)). 

Consider the L-operator 

^ ( t-:^ i^n^:z ) ■ 

It obeys the YBE (|^) if and only if A, C, D satisfy the commutation relations of the Uq{gl{2)) algebra ^ 
@, [|0|, |l|, H : 

AB = qBA, BD = qDB , 



DC = qCD, CA = qAC , 

A2 _ rj2 (4) 

q-q 

[A,D]=Q. 

This quadratic algebra has two central elements. One of them is a (/-analog of the Casimir operator: 

w = q-'^ A^ + qD^ + [q - q-^f BC , (5) 

another one, 

w^^AD, (6) 

for the Uq{sl{2)) case should be put equal to 1. If q is a root of unity some additional central elements 
appear. 

Irreducible finite-dimensional representations of dimension 2j + 1 can be expressed in the weight basis, 
where A and D are diagonal matrices: A — diag{q^ , . . . , q^-'). An integer or halfinteger j is spin of the 
representation. There exists the following realization by difference operators acting in the linear 
space of polynomials F{z) of degree 2j: 

AF{z) = q-^F{qz) , 

BF{z) = --^{q-^^F{qz) - q^'F{q-'z)) , 

(7) 

CFiz) = —AFiqz) - F{q-'z)) , 

z[q-q ^) 



DF{z) = q^F{q^'^ 



z 



Then Fq{z) = 1 is the lowest weight vector whereas F2j{z) ~ z^^ is the highest weight vector, i.e., 
CFq{z) = 0, BF2j{z) = 0. The Casimir operator (0) in this realization is equal to the c-number + 

If g is a root of unity there is, in addition, three-parametric family of finite-dimensional representations 
having, in general, no lowest and no highest weight Sometimes they are called cyclic representa- 

tions [^ . The difference quasi-exactly solvable equations corresponding to the cyclic representations of 
Uq{sl{2)) are particularly important [Q in applications to the Azbel-Hofstadter problem. However, in 
this paper we do not consider this case. 

2). L-operators associated to Aq{SL(2)). 

'^We use Koornwinder's notation ]24[ for the generators. 
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Another important class of L-operators is constructed by means of the dual quantum algebra Aq{SL(2)), 
the g-deformed algebra of functions on the group SL{2). Consider the operator matrix 



It is easily verified that it satisfies the YBE (|^), 

R{u/ v)gi {u)g2 (v) = 52 {v)gi {u)R{u/ v) , (9) 
if and only if a, 6, c, d obey the algebra 

db = qhd , bd = qdb , 
dc = qcd , cd = qdc , 

[d,d] ^ {q~q ^)bc , 

[b,c]=0. 

These are commutation relations for the generators of the dual algebra of Uq{gl{2)) |2^. We denote 
it Ag{GL{2)). The conventional interpretation of Aq{GL{2)) identifies it with a g-deformed algebra of 
functions on the group GL{2). 

The central element is bc^^ (it belongs to an extended algebra); another one is the g-determinant 
dd — qbc, which for the S'L(2)-case should be put equal to 1. The corresponding factoralgebra is de- 
noted Aq{SL{2)). Restricting to the compact real form of the quantum group, one obtains the algebra 
Aq{SU (2)) , which was extensively studied [^ . 

For completeness, we give the list of irreducible unitary representations of Aq{SU{2)) for real q ||2^ . 
There are two series: 

a) . One-dimensional representations: 

6 = c = 0, d = d-^^e"^, (11) 

<ip <2tt. 

b) . Infinite-dimensional representations ( par ametrized by the same "continuous spin" ip. They can 
be realized on functions in one variable ||2^ , [gs) : 

df{z)^-z-'{f{qz)~f{q-^z)), 
bf{z) = q-^e^^f{qz), 
cf{z) = -q^e-'"^f{qz) , 
df{z) = qzfiqz) . 

Note that e**^ enters (|lj) in the same way as the spectral parameter u enters (^). 
There are also some c- number solutions to (|l|): diagonal, 

and antidiagonal, 

^(-)-(* 0' ^''^ 

where * stands for an arbitrary u-independent c-number. The former come from the i-operators g{u) 
taken in the representation of the series a) ( [Til ) while the latter correspond to the L-operator (^) and the 
one-dimensional representation {A = D = 0, B and C are c- numbers) of the algebra (|4|). 



(12) 
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The trace (in the auxihary two-dimensional space) of T{u) obeying (|^) is a generating function of 
commuting integrals of motion: t(u) = TrT(u), 

[tiu),t{v)]=0. (15) 

In the case of elementary L-operators there is only one independent integral of motion (considered as 
a hamiltonian). Though the commutativity ( |l5|) is meaningless in this case, the transfer matrix t{u) 
possesses all necessary formal properties, which allow one (at least, in principle) to apply the technique 
of the algebraic (or functional) Bethe ansatz. 

The rational limit means q — e^\u = e''", ?i — ^> 0, and u becomes an additive spectral parameter. For 
future reference, let us present some formulas related to the rational limit. In the rest of this section we 
write simply u instead of u. 

The L-operator (||) becomes 

where Si are generators of sl{2): 



[S±,So]^TS±, [S+,SJ^^2So. (17) 

The correspondence with Uq{sl{2)) is as follows: {A - D)/{2h) ^ Sq, B ^ S+, C ^ S-. 

The realization (^ is a smooth g-deformation of the standard representation oof si (2) by first-order 
differential operators: 

S- = ^.. So = z^-j, S+ = -z^-^^+2jz. (18) 

This representation has been used to constract and classify linear differential equations having polynomial 
solutions 1^, [Q, |l^. The QISM interpretation is given in Section 5. 

The rational limit of (||) is simply a (w-independent) generic "group element" of SL(2) taken in the 
fundamental representation: 

9iu)^[l J). (19) 

Note that a, b, c, d are c-numbers in this case since the algebra (|l^) becomes commutative. In Section 
6 this "L-operator" is used for a QISM interpretation of global S'L(2)-rotations of quasi-exactly solvable 
hamiltonians. 



3 General properties of monodromy matrices for open inte- 
grable spin chains 

Here we give a brief summary of the formalism treating integrable systems with boundaries. The boundary 
conditions consistent with integrability are determined by c- number 2x2 matrices Ki{u) and Kr{u) 
(for the left and right boundary respectively) depending on the spectral parameter and satisfying the 
"reflection equations" (RE) fT^I, 

R{u/v){Ki{u) ® \)R{uvq-^){\ ® Ki{v) = 

(20) 

= (1 ® Ki{v))R{uvq-^){Ki{u) ® \)R{u/v) , 
R{v/u){Kl{u) ® l)R{u-^v-^q-^){l ® Kl{v)) = 

(21) 

(1 ® Kl{v))R{u-^v-^q-^){Kl{u) (g> l)R{v/u) 

{t means the transposition) with the i?-matrix (0). Solutions for Ki and are related: if Ki(u) is a 
solution to (pO|), then Kf{u~^) is a solution to (|2l|). In the scattering picture, the RE's describe the 
factorized scattering of a two-state particle on left and right walls respectively. 

More generally, one can consider operator solutions to the RE's (|20|), ( [2l| ) (i.e., matrices K{u) with 
operator valued matrix elements). Speaking informally, the wall may carry quantum numbers. In this 
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case the RE's determine commutation relations of the matrix elements. They generate the "reflection 
algebra" . In the case of absence of the spectral parameter this algebra was studied in |^ . 

The QISM approach to integrable systems with boundaries was developped by Sklyanin The 
main results of the paper JTsf are summarized below in the form of two theorems. 

Theorem 1. Let T{u) satisfy the YBE (|l]) with the i?-matrix ^ and let Ki{u) (resp., Kr{u)) satisfy 
the RE (resp., (|2^)) with the same i?-matrix. Then 

Ki{u) = T{u)Ki{u)a2T'{u-^)a2 , (22) 

lCr{u) = {T\u)Kl{u)a2T(u~^)a2)' . (23) 

satisfy (^) and (^ij) respectively (here and below ai are Pauli matrices). 

Remark 1. The theorem holds for both operator and c-number ii'-matrices. 

Remark 2. One may interpret (p^), ( p3| ) as a "dressing transformation": K "dressed" by T yields IC. 
Remark 3. For unimodular c-number matrices independent of u the operation CT2r*cr2 is simply T^^ . 
Remark 4- For a c-number matrix the (operator) matrix K, is called the quantum monodromy 
matrix for an integrable system with non-periodic boundary conditions. 
It is convenient to represent (p^), (p^ pictorially as follows: 



ICi{u) = 




T T 



Theorem 2. Let K.i{u) and K,r{u) be any solutions of (|2C|) and ( |2lD respectively. Then the quantities 

t{u) = Tr(/C.(M)/C/(w)) (24) 

form a commutative family: 

[t{u),t{v)]=Q. (25) 

The quantity t(u) is called a quantum transfer matrix. Its diagonalization can be performed by 
means of the algebraic (or functional) Bethe ansatz technique. To describe integrable open spin chains, 
one should put ICr — Kr (a c-number solution) in (M) and substitute /C/ from (22): 



t{u) = Ki 




Kr 



4 Trigonometric case 

The boundary matrices for the reflection of a two-state particle on a scalar wall are given by js^ 

^'^''^-\x^{q-^u'-qu-^) -2xoisu-s-^u-^) J ' 
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''^ '~\y-{qu^-q-^u-^) -2yo{t-^u-tu-^) ) ' ^^^^ 

where Xo, x±, s and yo, y±, t are arbitrary parameters characterizing the boundary conditions. 

The L-operator L{u) (|) still satisfies the YBE (0) if u is multiplied by a constant k. Substituting 
L{uk) for T{u) in ( p^ we get the quantum monodromy matrix 

M{u) = L{uk)Ki{u)a2L\u-^k)(j2 (28) 

(in this specific case we denote it M{u)). The calculation of its matrix elements is straightforward. One 
should take into account that the quadratic Casimir element w (^) under any irreducible representation 
acts as a c-number. It is convenient to represent the result in the following form. 
Consider the operators 

Hi = x+kAB + x-k~^CA + 2{q - q'^y^xos'^A^ , (29) 

H2 = x+k-^DB + x^kCD - 2{q - q-^^XosD'^ , (30) 

H3=iq - q-y^x+ik^A^ + k-^D^) - {q - q-^)x^C^ + 2xoisk-^DC - s'^AC), (31) 

SsHq - k'^D^) - (q - qyx+B^+2xoiskBD - s-^k-^BA). (32) 

We note that Hi , H2 , H3 form a simple quadratic algebra (a slightly different version of this algebra was 
previously studied in ||3l|| ): 

q^^HaHfj - qHfjHa = 57-^7 + ^7 > (33) 
where {a, /?, 7} stands for any cyclic permutation of {1, 2, 3}. The structure constants are: 

gi^ -{l + q-^)x+, g2 = -{l + q^)x+, 9z = [q + q-^)x_ , (34) 
hi = 2-^2^(s(fc2 + fc-2) + q-^s~^w) , 

h2^-2^^^^{s-\e + k-^)+qsw), (35) 

/13 = ^——rUxl - x+x^ik^ + k^^)w) . 

q-q 

The operators Hi, H2, H3 form a similar algebra; in particular, 

q-^H2Hi - qHiH2 = iq + q-^)x+H^ + /13 • (36) 
Let us decompose M{u) into the operator part M{u) and the c-number part M^^^(m): 

M{u) = M{u) + A'&^ {u) . (37) 



Then M{u) can be compactly written down in terms of the operators (29)-(p2[): 

.y, . q^u^-qu-^ f -uHi + u-^H2 H3 \ 

^'^""^^ q-q-^ \ H, q-^U-^Hi-quH2) ^''^ 

Note that the boundary parameters do not appear explicitly in ( |3^ ) entering only through the structure 
constants of the algebra (BSh . For the c-number part one has 



M[l\u)^A4fiqu-^)= ^^ Msu-s-^u-^)Hk-^ + k-^^^^^^^^ 

(39) 
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It is clear from (|3l|), ( ^2|) that the non-diagonal elements of Af (w) generally do not have a zero mode 
(a "false vacuum") independent of u. In such a case the general strategy of the functional Bethe ansatz 
consists essentially in passing to the new basis formed by the eigenvectors of or . Note that these 
operators contain only elements of the lower (resp., upper) Borel subalgebra of Uq{sl{2)) (i.e, for example, 
i?3 is a quadratic form in A, D and C, not B). This allows one to find the eigenvectors of and H^, in a 
quite explicit form Q . Under the representation the eigenfunctions are big g- Jacobi polynomials (see 
e.g. ^^). This fact may be useful for diagonalization of t{u) by means of the functional Bethe ansatz. 
Disregarding the c- number part, we get the quantum transfer matrix: 



t(u) = Tr{Kr{u)M{u)) = 



-{2yo{t-'H2-tH,)+y+H3+y_H3). (40) 



It is clear that in this case the family of commuting integrals of motion generated by t{u) contains only one 
(independent) operator. In terms of A, B, C, D (^ the transfer matrix becomes a generic homogeneous 
quadratic form in these operators (see (p9|)-(p^)) since it depends on 7 parameters: 3 in each iiT-matrix 
(the common factor in (^6|) is inessential) and k. Indeed, the total number of coefficients of a general 
quadratic form is 10 but two of them contribute only to the c-number term in ( p7| ) due to the two central 
elements (AD — 1 and the Casimir operator); besides, the common multiplier is also inessential. 

As it is shown in the hamiltonian 11 of the Bloch particle in a magnetic field is a particular 

quadratic form in the Uq{sl{2)) generators with \q\ = 1 (the coefficients depend on the gauge and the 
type of the lattice) and therefore this system can be considered as an integrable model. Here is a list of 
the most important examples [0, 1^, [Q. 

1) . Square lattice, modified Landau gauge: Ki = CT3, Kr = cri, fc ^ 00, 

H = -i{q~q-^)q-^/^{CA + BD). (41) 

The flux per plaquette is $ = 2-!tP/Q (P, Q are coprime integers), q = e'*/^. 

2) . Square lattice, chiral gauge: Ki = + {q~^^^u — q^^^u~-^)a^, Ki ^ + {q^^^u — q^^^-^u^^)a-, 
k = l, 

H = i{q- q-^)q-^/^{CA - BD + qBA - qCD). (42) 

The flux per plaquette is AttP/Q [P odd). 

3) . Triangular lattice, modified Landau gauge: ^ ~ 0, xq — yo = k — 1, t = —s ^ 00, 
X- ~ y^ 00, y+/s —>■ 2Ag^/^ exp(— i7r(P — l)/2), where A is the hopping amplitude along the third 
axis, the flux per elementary triangle is ttP/Q {P odd). The hamiltonian is 

H = A(yl2 + D'^) + e''^(-P-i)/2g-i/2(q _ q-^){CA + BD). (43) 

4) . Square lattice, chiral gauge (a different version). We include this example only for completeness. It 
is related to transfer matrix for a closed chain: t{u) = Tr{a2L{u)) (see ([T^). Here L{u) is the L-operator 
(^) and a2 is a c-number solution to (|l|) of the form (|lj). The hamiltonian is 

H^l{q-q-'){C-B), (44) 

{Q odd, P even). 

Under the representation (|^) the transfer matrix ( ^0| ) becomes a second-order difference operator in 
z. Clearly, this operator has the invariant subspace of polynomials spanned by 1, z, z^, . . . ,2:^-'. The 
polynomial eigenfunctions lying in this "algebraic" sector and the eigenvalues can be found in the form 
standard for the algebraic Bethe ansatz technique. A detailed analysis of the equations arising from ( ^0| ) 
and their Bethe ansatz solutions is given in Q (see also Appendix B in [D). 

Here is the explicit form of these equations in terms of the parameters of Ki and Kr- By means of 
(0) and (|2|)-(|2|) we rewrite the spectral problem for (|o|) as follows: 

a{z)^p{q'^ z) + d{z)ip{q^'^z) - v{z)'iIj{z) = Etpiz) , (45) 

where 

a{z)=(q-'^^+'^x+z~^q'^ ^x^k-^z-^(-q-^hj+z + 2yotkq-^ +q-^y^k^z-^^ , (46) 
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v{z) — q^^a{q' z) + q^'^^ d{q^^ z) — 



yo2:-(— ] — ) + y-Xa{ -, ) z 



(48) 



t 

Note that the dependence on the left and right boundary parameters in a{z) and d{z) completely factor- 
izes. 

In general, some eigenfunctions of ( ^ ) lie beyond the algebraic sector. Hence this spectral problem is 
quasi-exactly solvable, i.e., only a part of the spectrum is available in a closed algebraic form. There is one 
very important exception, where the spectrum itself is finite and the algebraic sector totally covers it. This 
is the case when g is a root of unity {q — e<-^PIQ ^ as before) and j = {Q — l)/2, so we obtain an important 
class of periodic difference equations on a finite ring-like lattice. Then all periodic eigenfunctions can be 
found algebraically Q. Furthermore, all quasiperiodic eigenfunctions (corresponding to generic internal 
points of bands) are expressed in terms of the family of the cyclic representations of Uq{sl{2)). This is 
just the case relevant to the Azbel-Hofstadter problem. In particular, for Example 3) above one obtains 
the spectral equation 

(z'l - \q)tl){q^z) + {q-'^z - \q-^)i){q-'^ z) - (z + z-^)tl){z) = E^l^{z) . (49) 

which provides the midband points of the spectrum Another approach to the Azbel-Hofstadter 
problem on a triangular lattice was proposed in p4||. 

Let us point out an alternative way to convert ([40| ) to a difference operator. One may disregard the 
explicit formulas (|2^)-(32|) connecting Hi with A, B, D and make use of the representation theory 



of the algebra (|33|), (p6) developped in Namely, the generators Hi, H2, H^, H^ can be realized as 
second-order difference operators of an essentially different form than the one following from (p9|)-(|3^ 
and (^ . It is evident from (^) that under proper conditions the two forms are equivalent (at least in the 
algebraic sector). More precisely, the corresponding operators, being restricted to the algebraic sector, 
are connected by conjugation with a certain matrix, i.e., by an isospectral transformation. In some cases 
such a transformation exists for the difference operator itself, not only for its algebraic truncation. In this 
paper we will not discuss this interesting question and only mention the following specific example. The 
eigenfunctions of H^ ( pl|) under the representation (0) (where A and D are diagonal) are big g-Jacobi 
polynomials (see e.g. p2[). On the other hand, as it follows from ||3l|, the eigenfunctions of H^ in the 
basis, where Hi is diagonal are Askey- Wilson polynomials ||3^. Our results indicate that there should 
exist a similarity transformation between the two difference operators. Some more details may be found 
in |. 



5 Rational limit 

In this section we consider the limit 9 — > 1 (the "rational", or continuum limit) of the monodromy matrix 
(^), providing a basis for embedding the continuous quasi-exactly solvable problems ||l^, ||l5| into the 
quantum inverse scattering approach. 

The construction of the previous section may have different continuum limits. Below we consider the 
most important one, which is directly related to integrable models with the rational i?-matrix. The rule 
of performing this limit is as follows. Put q — e/^, u = e''", s — e''*, i = e'**, k = e^^, and find the h- 
expansion of the monodromy matrix ( p8|) as ?i ^ 0, provided xq, x±, yo, y± are ?i- independent constants. 
In doing so, we will write u, s, t, k instead of it, s, t, k respectively. Since the trigonometric and rational 
cases never mix, this convention can not lead to a confusion. Note that u becomes an additive spectral 
parameter. 
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This rule is equivalent to repeating the general arguments of Sects. 2 and 3 in the context of the 
rational i?-niatrix and the L-operator (|l^). The general c- number solutions to the rational RE are |^ 

/ 2xo{u-s-l) x+i2u-l) \ , . 

'"^^ x^{2u-l) -2xoiu + s) J' ^ ' 

( 2yo{u + t + l) y+{2u+l) \ 
'■"^ y-(2.i + l) -2yn{u-t))- ^"'^ 

However, it is more convenient to proceed by considering the limits of (|^ 
The operators (p9|)-(^2|) are expanded as 

Hi = n-^xo + h["'> + hH[^'> + o{n^) , 



H3 = h'^x+ + 2hH^^'' + o{n'^) , 



(52) 



where the operator coefficients are expressed through the gener ators (|l^) of sl{2): 

h[°'> = i/f ^ = 2a;oS'o + x+S+ + x^S- - xqs , (53) 

h[^'> = -i?^ ^ = 2xoSf) + x+SqS+ + X-S-So - 

-2xqsSq + x+kS+ - X-kS- + -^xais^ - -) , (54) 

H^^^ = x+Sl - x-Sl - 2xoSoS- + 2x+kSo + 2xo(s - k)S- + x+{k^ - ^), (55) 



h''^^ ^x^Sl- x+Sl - 2xoS+Sq - 2x^kSo + 2xq{s + k)S+ + x^{k^ - . (56) 



1 

12' 

The operator M{u) ( ^8| ) aquires a c-number part as ?i ^ 0, with the leading term being singular (~ ?i~^). 
In what follows we neglect all next-to-leading c-number contributions since they are absolutely irrelevant. 
In particular, we can throw away the c-number terms in (p3|)-(p6|). Moreover, it is easy to see that the 
?i~^-terms exactly cancel in the sum (|37|). Finally, one obtains the following rational monodromy matrix: 



^^""^H ifW if«-(.+l)i/(") 



(57) 



Combining it with (|5l|), we get the transfer matrix: 

t{u) = (4«2 _ 1) (y+J?(i) + y_Hf^ - 2yoi/f ^ - 2yoti?}"^) = (58) 

= (4u^ - 1) {{y+x,- + y-xj^ - 4?;o2;o)'S'o - V+x^S\ - y-x^St- 

-{y+xo + yox+){S+SQ + SqS+) - {y-Xo + yoX-){S-SQ + SqS-) + 
+2{k{y^x+ - y+x^) + 2{s - t)yoXo)So + 
+2{{s + k + l/2)y+xo -{t + k+ l/2)yox+)S+ + 

+2{{s - fc + l/2)y_xo -{t-k + l/2)yoX-)S- ) . (59) 

It is a generic mixed quadratic- linear form in Si. The number of independent parameters is the same as 
in (Eol). 
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The diagonalization of (|5^) gives (for the spin j representation the following differential equation: 



dz^ V 2'^ ^ ' '7 dz 



(60) 



where 

Q{z) = {x+z"^ - 2xqz - a;_)(y+z^ - 2yoz - ?/_), (61) 

P{z) = -{{s + k + l/2}y+XQ- {t + k + l/2)yQX+)z^ + 
+ {k{y-x+ - y+x-) + 2{s - t)yoXo) z + 

+ {s-k + l/2)y^xa - {t - k + l/2)yox^ . (62) 



Equations of this type are well-studied. If Q{z) has 4 simple roots, the eigenvalue problem (60) can be 
reduced to Heun's equation [3^ . The transformation of ( |60| ) to the Schrodinger form is discussed in detail 
in the reviews ]l3|, |l4[ 



A relation of (30|) to integrable spin chains was already mentioned in the literature |15). However, the 
known relation is absolutely different: the operator in the l.h.s. of (^0|) is identified with the hamiltonian 
of inhomogeneous Gaudin's magnet psf on 3 (or 4) sites with periodic boundary conditions. It is known 
that Gaudin's magnet is a quasiclassical limit of the integrable spin chain with the rational i?-matrix. 

In the present paper we identify ( |60|) with the eigenvalue problem for the transfer matrix of a formal 
XXX-typ^ "magnet" on only one siten but with non-periodic boundary conditions. As it is shown in 
Sect. 4, the trigonometric generalization of this system leads to a class of quasi-exactly solvable difference 
equations whereas an analog of the Gaudin's magnet picture for the latter is not known. 



6 Adjoint action of SL{2) 



There is an obvious group of isospectral transformations of (58). These transformations are induced by 
the adjoint action of SL{2): Si g~^Sig. One may say that really different spectral problems correspond 
to S'L(2)-orbits in the space of quadratic forms (^). The generic orbit is 3-dimensional, so the number 
of independent parameters is reduced to 4. Under the adjoint action of a group element 



a b 
c d 



ad — be = 1 , (63) 



the generators transform as follows: 

S+ ^d^S+ + 2edSa -c^S-, 

So bdS+ + (1 + 2^)5*0 - aeS- , (64) 

S- -b^S+ - 2abSQ + a^S- . 

Making this transformation in (|58|), one gets an operator having the same spectrum for any spin j. Our 
aim in this section is to show how this transformation may be interpreted in terms of the QISM. 

To do this, recall the c-number solutions ( p^ of the rational YBE. Let us consider the group element 
( |63| ) as such a solution and apply Theorem 1 (see Sect. 3) to the pair g, Ki (or g, Kr) for Ki (Kr) given 
by ( |50| ) ((|5T|)). It is convenient to fix common multipliers in ( |50| ) and (^ij) by putting (2s + l)xo — 
1, {2t + 1)2/0 — 1- Then Ki (resp., Kr) depends on a 3-dimensional vector x = (xi, X2, xq) (resp., 
y — (2/1, 2/2, yo))i where x± = xi ± X2, y± = yi ± y2 ■ Working in this normalization, we indicate the 
dependence on x and y explicitly: 

Ki{u; x) = {2u - 1) (xa) - 1 , (65) 

*This "spin chain" has only one spin because ( p8[ ) includes just one elementary L-operator. 
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Kr{u;y)^{2u+l)iya) + l, (66) 

where (xtr) = xi<ti + X2<72 + xoa^ denotes the inner product of 3-dimensional vectors (here cr^ are PauU 
matrices), and 1 is the unit matrix. 

Now, according to Theorem 1, we should "dress" if-matrices using (|2|), (|3|). In the simple case at 
hand the dressing is reduced to the conjugation: 





where the dashed lines (carrying trivial one-dimensional "quantum space") denote the insertions of the 
"L-operator" g. The equalities immediately follow from (65), (pq). The shorthand means the adgoint 
action of g to the 3-component vector x. One concludes from (|67|) , ( |6^ ) that the dressing in the case at 
hand is equivalent to a rotation of the vector parameter. 

Let us represent the generators Si as a vector with operator-valued components: 



ihs^ + s+), hs^-s+), So). 

2 Zi 



(69) 



The adjoint action g is given by (p^. In these terms the transfer matrix (58) can be written in the 
form: 



t(u) = (4u2-l)(((yxS)(xxS)) + ((xxS)(yxS))-(yS)(xS)- 

- (xS)(yS) + 2((y - x)S) + Aik{{y x x)S)) + c-number , (70) 

where a x b denotes the skew product of 3-vectors: (a x b)^ — eaf}-yO,pb-y . It is clear from ( |70| ) that the 
operator part of t(u) is invariant under simultaneous rotations of all the vectors x, y and S. In other 
words, the rotation g^^S ( |64|) in ( [70|) is equivalent to x — > gx, y ^ gy given by the "dressing" (|67|), (|6^). 
The dressing means the insertion of g to the left and g^^ to the right of the line corresponding to L{u) 
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(|l6|). Schematically, 




Another way to see this is to observe that 

g-^L{u + k)g ={u + k)l + {(j{g-^^)) (72) 

for the L-operator (|l6|), and (j2L*(—u)a2 = —L{u) . Then 

t{u) = -Tr {Kr{u] y)L{u + k)Ki{u; x.)L{u - k)) , (73) 

and the transformation S .9^^S in L{u ± k) leads to 

T^'\u) = -Tr{Kr{u-y)g-^L{u + k)gKi{u-^)g-^L{u~k)g) = 

= -Tr {Kr{u;gy)L{u + k)Ki{u;gx)L{u-k)), (74) 

which is equivalent to ( |70| ) due to (|67|), (|68|). 

It is an interesting open problem to find a proper q-analog of the transformation considered in this 
section. In particular, it is not known whether there are any isospectral subfamilies among the operators 
of the form ( ]40| ) other than the trivial ones {B — > e~-^'"^B, C e^'^'^C), which correspond to the similar 
insertion of g{u) (||) taken in the one-dimensional representation (|ll|). We hope that our approach may 
help to solve this problem. 



7 Concluding remarks 

We have shown that both discrete and continuous quasi-exactly solvable problems of quantum mechanics 
are tractable in the framework of the quantum inverse scattering method. Quantum transfer matrices 
for a peculiar simple integrable system with boundaries yield the complete collection of quasi-exactly 
solvable hamiltonians. These hamiltonians are quadratic forms in the generators of Uq{sl{2)) (or sl{2) in 
the rational limit) taken in a finite-dimensional representation. 

This reformulation opens a way to apply the powerful methods specific for quantum integrable systems. 
For the representations of Uq{sl{2)) having both highest and lowest weights these methods give the results 
which are eigther already known or can be obtained by means of more elementary tools Q. However, 
if g is a root of unity there exists a family of cyclic representations having in general neigther highest 
nor lowest weights. This is just the case relevant to the Azbel-Hofstadter problem, where the generic 
points of bands are described by cyclic representations. The Bethe ansatz in this case is much harder 
problem. Here the reformulation in terms of the QISM may led to really new outcomes. The appropriate 
method is the technique of Baxter's intertwining vectors applied to the chiral Potts model (which is also 
known to be connected with cyclic representations) in Q. Recently, this method was applied Q to 
the Azbel-Hofstadter problem. Presumably, the method should work for any operator of the form ( ^0| ) 
as well. 

Let us point out two questions, where the results of the present paper may contribute something to 
the conceptual understanding. 

One of them was already mentioned at the end of Sect. 6. It is the question about isospectral subfamilies 
among operators of the form ( ^ ) (or, rather, about a proper g-analog of them). In the continuous case, 
they are orbits of the adjoint action of SL{2) group elements in the space of quadratic forms in sl{2)- 
generators. We have seen that the group elements may be considered as L-operators (obeying the YBE 
with spectral parameter), the adjoint action being an insertion of them into the monodromy matrix. 
Remarkably, each ingredient of this picture has its natural g-deformed counterpart. The notion of a 
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quantum group-like element was recently discussed [ p7[ (from another point of view) in connection with 
quantum r-functions. A comparison of these studies with our results may be fruitful for both approaches. 

A related question concerns separation of variables. It is known |Q that inhomogencous n-site 
Gaudin's magnets are in one-to-one correspondence with separated coordinate systems for the Laplace- 
Beltrami operator on the {n — l)-dimensional sphere (or hyperboloid). On the other hand, Gaudin's 
magnet on 3 sites generates continuous quasi-exactly solvable hamiltonians [ p^ . Considered as quadratic 
forms in generators of sl{2) such a hamiltonian determines a separated coordinate system on the 2- 
sphere (or hyperboloid). The non-equivalent separated systems correspond Q to S'L(2)-orbits (under 
the adjoint action). The quadratic forms in generators of Uq{sl{2)) might have a similar relation to 
hypothetical "separated coordinate systems" on quantum spheres and hyperboloids. 
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